function output =  testODEfunc(t, y)
    arg = 2;
    dy = zeros(3,1);
    
    dy(1) = arg*y(1) + y(2);
    dy(2) = -arg*y(2)+2*sin(y(1));
    dy(3) = y(3)-y(2);
    output  = dy;
    
    
    return